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Determining the upper limit of Tee for the 1^(4260) 

X.H. Mo G. Li C.Z. Yuan K.L. He H.M. Hu J.H. Hu P. Wang Z.Y. Wang 
^Institute of High Energy Physics, P.O.Box 918, Beijing 100049, China 
''China Center of Advanced Science and Technology, Beijing 100080, China 

"^College of Physics and Information Technology, Guangxi Normal University, Guilin 541004, China 

By fitting the R values between 3.7 and 5.0 GeV measured by the BES collaboration, the upper limit of the 
electron width of the newly discovered resonance y (4260) is determined to be 580 eV/c^ at 90% C.L. Together with 
the BaBar measurement on ree B(y(4260) ^ tt+tt^ J/t/;), this implies a large decay width of ^(4260) n^n^J/tp 
final states. 



1. Introduction 

Recently, in studying the initial state radia- 
tion events, e+e~ — > jjsrtt^t^^ J/ip {jisR- ini- 
tial state radiation photon) with 233 fb~^ data 
collected around ^/s = 10.58 GeV, the BaBar 
Collaboration observed an accumulation of events 
near 4.26 GeV/c^ in the invariant-mass spectrum 
of TT~^Tr~J/tjj The fit to the mass distribu- 
tion yields 125 ± 23 events with a mass of 4259 ± 
8^6 MeV/c^ and a width of 88 ± 23^^ MeV/c^. 
In addition, the following product is calculated 

r(r(4260) ^ e+e-) • 6(r(4260) ^ tt+tt- J/V')) 

= 5.5±l.0ta.7 eV/c2. 

Since the resonance is produced in initial state 
radiation from e+e^ collision, its quantum num- 
ber J^^ — 1 . However, this new resonance 
seems rather different from the known charmo- 
nium states with J^^ = 1 in the same mass 
scale, such as i/;(4040), ?/;_(4160), and V(4415). 
Being well above the DD threshold, instead 
of decaying predominantly into final 
states, the y(4260) shows strong coupling to the 
TT^TT^J/ip final state. So this new resonance does 
not seem to be a usual charmonium state. The 
strange properties exhibited by the y(4260) have 
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triggered many theoretical discussions P]-|12|. 

One suggestion is that the y(4260) is the first 
orbital excitation of a diquark-antidiquark state 
([cs][cs]) |2I3| . By virtue of this scheme, the mass 
of such a state is estimated to be 4.28 GeV/c^, 
which is in good agreement with the observation. 
A crucial prediction of the scheme is that the 
F(4260) decays predominantly into DgDs- 

Another opinion favors a hybrid explanation 0J 
IFiGii?) . In the light of the lattice inspired fiux- 
tube model, the calculation shows that the de- 
cays of hybrid meson are suppressed to pairs 
of ground state 15' conventional mesons |13I14| . 
This implies that decays of r(4260) into DD, 
DgDs, and D*D* are suppressed whereas D*D 
and D*Ds are small, and D**D, if above thresh- 
old, would dominate (P-wave charmonia are de- 
noted by D**). So it is interesting to search for 
the possible decay of ^(4260) Di (2420)5. 

The third interpretation we wish to mention 
is provided by Ref. which suggests that the 
F(4260) is the second most massive state in the 
charmonium family. The author ascribes the lack 
of F(4260) in e~^e~ hadrons to the interference 
of S-D waves, and also estimates 



r(r(4260) 



) ~ 0.2 - 0.35 keV/c2 



(2) 



Besides the above interpretations, there are 
other kinds of proposals. The lattice study in 
Ref. 9 suggests that the y(4260) behaves like a 
DiD molecule. In Ref. [TOj, it is proposed that 
the new state might be a baryonium, containing 
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charms, configured as A^-Ac. In Rcf. TP, the 
y(4260) is considered as a p-Xci molecule while 
in Ref. the F(4260) is considered as an uj-Xd 
molecule. However, all aforementioned specula- 
tions need further experimental judgment. 

Most recently, CLEOc collected 13.2 pb^^ data 
at ^/s = 4.26 GeV and investigated 16 decay 
modes with charmonium or light hadrons |15| . 
and the channels with more than 3a statisti- 
cal significance are tt+tt^J/?/; (He), Tr'^ir'^J/ip 
(5.1cr), and K+R-J/ip {3.7a). No compelling 
evidence is found for any other decay modes for 
the r(4260), nor for the ^'(4040) and V'(4160) res- 
onances [15] . These measurements disfavor the p- 
Xci molecular model ^J, baryonium model |1U| . 
and high charmonium state explanation So 
far as other surviving speculations are concerned, 
such as charmonium hybrid |4I5I6I7| . tetraqark 
model |2I3| . DiD molecule suggestion, and uJ-Xci 
molecule explanation }1 2j , further experimental 
studies are needed to make more definitive con- 
clusions. 

Since the y(4260) was observed in e+e^ anni- 
hilation, it is expected that it contributes to the 
total hadronic cross section in e'^e~ annihilation 
(or the R value, in other words). The most re- 
cent such data on R measurements are from the 
BES experiment jl8ll9| . as shown in Fig. ^ for 
Ec.m.{= Vs) from 3.7 to 5.0 GeV. If we look in 
detail within the range from 4.25 to 4.30 GeV (re- 
fer to the inset of Fig. it seems there is a bump 
around 4.27 GeV. Has this structure a connection 
with the y(4260)? We try to answer this ques- 
tion in this Letter. As there are other resonances 
nearby, we shall fit the full spectrum between 3.7 
to 5.0 GeV in order to get the information on the 
y(4260). 

2. Fit to the ^-resonances 

The R values |18I19| used in this analysis were 
measured with the Beijing Spectrometer (BESII), 
which is a conventional solenoidal detector ex- 
pounded in Ref. j^. In the analysis below, the R 
values are converted into a cross section a{e~^e~) 
by multiplying the Born order cross section of 
e+e" P-^ P~ ■ These cross sections are plot- 
ted in Fig. 121 There are clear peaks of i/'(3770). 




Figure 1. R values for Ec.m. from 3.7 to 5.0 GeV 
measured by the BES collaboration |lj^ll9) . The 
inset shows the R values in the vicinity of the 
r(4260). 



■0(4040), -0(4160), and 0(4415). The data points 
have been used to obtained the parameters of the 
■0- resonances [T7) . 

To acquire the resonance parameters, we could 
simply fit the data set with cross section formula, 
each resonance with a Breit-Wigner 



^ [S- (AP)2]2 + (MJT^')2 ' ' 

where Fee, F^, and Tt are the mass indepen- 
dent electronic, hadronic and total widths, re- 
spectively, for a vector resonance of mass M pro- 
duced in the head-on collision of e+ and e~. No- 
tice that Fee ^ r* in our analysis; the approxi- 
mation F/i w Ff is actually adopted hereinafter. 
The summation of all the cross sections within 
the range we are studying is 

4 

^^^]{s) , (4) 
where indices 1, 2, 3, and 4 denote four resonances 
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Figure 2. Total hadronic cross section in nb ob- 
tained as a{e'^e~ — > hadrons) = R ■ 86.85/s (s 
in GeV^) from R values in Rcfs. I 18ll9j . Three 
sets of fit results, Fit 1, Fit 2, and Fit 3 corre- 
spond to three cross section forms ai, a2, and 
CTa as described in text. Two interference curves 
have been moved downward by 5 nb for display 
purposes; the dashed line at —5 nb corresponds 
to zero cross section in the fit. 



^-(3770), V(4040), ?A(4160), and ^(4415), respec- 
tively. 

In the -i/j-family resonance region, if we assume 
that all two-body states are decay prod- 

ucts of resonance, and not produced directly in 
continuum, we could therefore treat resonance 
and continuum incoherently. Nevertheless, for 
the four wide resonances, they are close and are 
expected to have some same decay final states, 
there must be interference between any two of 
the resonances. Therefore the amplitudes corre- 
sponding to each resonance, with the following 



form 



i2^r^,rL 



(AP)2] +iAfJT: 



(5) 



have to be added coherently to give the total am- 
plitude that, once squared, will contain interfer- 
ences of the type ^T*Tj. If the resonances are 
quite broad, the interference effect can also dis- 
tort the resonance shape, the width might appear 
broader or narrower, and the position of the peak 
can be displaced as well. In this case, the total 
cross section is 



CT2(s) 



(6) 



where Tj (s) is given in Eq. (O . 

So far as the amplitude is concerned, in prin- 
ciple, there are presumably relative phases be- 
tween different amplitudes besides the phase due 
to complex Breit-Wigner formula itself. So a 
more comprehensive total cross section would be 
the summation of amplitudes together with an 
additional phase, viz. 



(7) 



Since what we actually obtain is the squared 
modulus of amplitudes, only three relative phases 
could be detected in practice. 

The standard chi-square estimator is con- 
structed as follows 



E 



^tht 



(Acre^P(Sj))2 



(8) 



where cr'^^P(sj) indicates the experimentally mea- 
sured cross section at the j-th energy point, while 
a^^'^isj) is the corresponding theoretical expecta- 
tion at this energy point, which is composed of 
two parts 



the 



is,) 



(9) 



where cr^°" denotes the contribution from contin- 
uum. Since there is little evidence in the data for 
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any substantial variation of the continuum back- 
ground within the studied energy region, we pa- 
rameterize the continuum cross section with a hn- 
ear function 



CT6g(s) = A + B{^s - 3.700), (s in GeV^ 



(10) 



as has been used in Ref. 17^ . Here we consider 
the continuum contribution as the background 
for measurement of resonance parameters, that 
is ab<,(s) =a™"(s). 

In Eq. © a'''^^ denotes the contribution from 
resonances. The fit resuhs are displayed in Fig.|21 
where Fit 1, Fit 2, and Fit 3 correspond to the 
three cross section forms ai, CT2, and 173, as ex- 
pressed in Eqs. 10} , ©, and iQ, respectively. 
Although the synthetic curves for three fits are 
almost the same, as we expected, the interfer- 
ence effect deforms each resonance significantly, 
according to Fig. |21 -0(3770) and V(4040) be- 
come narrower while V'(4160) and -0(4415) be- 
come wider when interference effects are included. 
There exist constructive interferences as well as 
destructive ones, and the interference behaviors 
for amplitudes with and without extra phases are 
also very distinct. The fit results indicate that 
Fee is very sensitive to the fit strategy, and the 
largest difference between various amalgamation 
strategies could reach 50%. 

Since 0(3770), 0(4040), 0(4160), and 0;(4415) 
have similar decay features, that is, all decay 
dominantly to final states, we prefer the 

synthetic scheme of amplitude with phase which 
takes into account all possible interactions be- 
tween resonances. 

So far as the y(4260) is concerned, the study of 
the BaBar Collaboration ^ implies that this new 
state may not be a common charmonium state, 
whose decay feature is rather distinctive from 
other ■0-family members, and this also obtains 
support from recent CLEOc measurements |15| . 
So it is favorable to perform the incoherent ad- 
dition for the contribution of the F(4260) to the 
total cross sections in Eqs. Q), or 

In addition, two other effects should be taken 
into account. First, there are theoretical ar- 
guments in some references |3l4l5lfi| that the 
y(4260) may be due to a threshold effect just as 
the 0(3770), so the width of the y(4260) could 



be energy dependent as follows: 

r_D_D(s) =Tdd ■ d{y/s ~ 2m^,±)y. 

1 + (rpjj.±f 



(11) 



1 + {'rpjj.±Y 



T)3 



Here the subscript "DD" explicitly denotes the 
threshold effect, Fdd = TDDiM"^), r is the clas- 
sical interaction radius, p is the D*^ momentum, 



Pd' 



s — Ami 



d: 



and p is the D*^ momentum at resonance peak, 
viz. 



1 



M2 - 4m2 .± 



Second, the F(4260) has been observed decay- 
ing into the 7r+7r~J/0 final state and this kind 
of decay can be expressed by the amplitude with 
the form of Eq. ij^l which is different from that 
of the decay into D*^D*^ with threshold effect. 
Furthermore, although the present data display 
the distinctive feature of the y(4260) from other 
charmonium resonances, some connections pre- 
sumably exist which could lead to the interference 
effects of the F(4260) with other -0-family mem- 
bers. With all these considerations in mind, we 
split the amplitude of the F(4260) into two parts. 
One is the amplitude of the threshold part, which 
is defined as follows: 



Tdd{s) 



ee 

{s - M^) + iMTtis) ' 



(12) 



where the total decay width is defined as the sum- 
mation of two parts : 



DD 



(^) 



Here the first term is the energy dependent par- 
tial decay width for the threshold part while the 
second term is the energy independent partial de- 
cay width for the non-threshold part (denoted by 
FnoTi)- We introduce a factor / to indicate the 
ratio of T]j]j{s) to ^t{s) at the resonance peak, 
that is _ 



f = 



Ft 
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Notice the energy independence of r„o„ , we have 
fnon = (1 - /)rt, then 

rt{.s)=TDD{s) + il- .f)Tt . 

In this case, the amplitude of the non-threshold 
part is given by the following expression 



Tis) 



12^(1 - f)TtT,e 
{s - AP) + iMTtis) 



(13) 



Without the knowledge of /, we first set / — 0.5 
in the study in Sect. |21 and leave the variation 
effect of / to Sect.H 

By virtue of the above discussion, the synthetic 
cross section within the region studied takes the 
following form 
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■T{s)e 



+ \TnD{sf -(14) 



Under such a scheme, the upper limit of the pro- 
duction of the F(4260) in e^e~ annihilation is 
determined, as expounded in the following sec- 
tion. 

3. Determination of T^e of the y(4260) 

Various fits to the data tell us that with 
limited knowledge on the nature of the reso- 
nances and comparatively meager data, we could 
only determine the T^e of the F(4260) by a 
scan method. Specifically, we fix the mass of 
the y(4260) at 4.259 GeV/c^ measured by the 
BaBar Collaboration and scan over the Tt G 
(20,180) MeV/c2 (8 MeV/c^ step) and e 
(0, 1000) eV/c^ (10 eV/c^ step) parameter space. 
In order to avoid some grotesque fit results due 
to random effect of the F(4260) on the reso- 
nances nearby, in the scan procedure, the lower 
and upper bounds of the masses (total widths) 
are fixed to be 4100 and 4220 MeV/c^ (30 and 
250 MeV/c2) and 4350 and 4500 MeV/c^ (30 and 
300 MeV/c^) for ^(4160) and -0(4415), respec- 
tively. 

For each pair of T\^{= [i x 10] eV/c^) and 
r|(= [20 + i X 8] MeV/c^), fitting the R data 
with the determined from Eq. ||SJ), we obtain 



a best estimated as a function of the r*„ and 

- t : 

1 

r 



X 

r^, or equivalently a relative likelihood, viz. 



£.(rL,n) = exp(-^x'(rL,n) 



(15) 



Instead of using this directly, we further con- 
struct a weighted likelihood as follows 



= /Ar-C(rL,n)x 

1 / {n 



■ cxp 



(16) 



where /jv is an arbitrary normalization factor 
and the Gaussian term indicates that the pos- 
sible total width (F^) of the r(4260) is consid- 
ered to distribute as a Gaussian with the mean 
value Tt — 88 MeV/c^ and the standard deviation 
= 23.8 MeV/c^, which have been determined 
by the BaBar Collaboration ^J. The weighted 
likelihood (Cw) as a function of the T^^ and 
of the r(4260) is shown in Fig. I^a). bumming 
£„(Fgg,F() with respect to F^, we obtain the 
variation of versus Fee as shown in Fig.OJb). 
The integral of the likelihood curve gives the up- 
per hmit of Fee of the r(4260) at 90% confidence 
level (C.L.): 



(17) 



By virtue of Fig. I^fa), we sum up £u,(Fee,Fj) 
with respect to Fge to 90% fraction of the total 
area then obtain the upper limit of F^^^^^"-* at 
90% C.L. for each rf^'^^^°\ So we obtain the 
variation of the upper limit of Fee versus F( as 
shown in Fig. 0] All the points are almost in a 
straight line, which indicates the ratio of the two 
quantities, or the upper limit of the branching 
fraction of y(4260) e+e^ does not depend on 
the total width. Taking the slope of the solid line 
in Fig. 01 we get 



6(y(4260) 



) < 4.6 X 10 



-6 



at 90% C.L. 



4. Other possibilities 



Although it is reasonable to treat all the ip- 
resonances coherently as we once mentioned in 
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(a) Lw w.r.t. Fee v.s. Ft 




Fee (eV/c') 

(b) w.r.t. Fee 

Figure 3. The weighted hkehhood distribu- 
tion with respect to (w.r.t.) Fee and Ft (a); and 
to Fee with Ft integrated (b). The "N.L" in (b) 
indicates normalized integral value for L^. 



Sect. 121 one may argue that other possibili- 
ties could exist since there is limited informa- 
tion about the properties of these ^/^-resonances. 
Thereby it's better to take a variety of effects into 
consideration. 

First, there are three schemes for summation 
of '0 resonances in this analysis: through cross 
sections, amplitudes, or amplitudes together with 
relative phases, as discussed in Sect. |21 the cor- 
responding results are Fee < 30, 10, 420 eV/c^ 
at 90% C.L. respectively, where the last one has 
been given in the previous section. 

Second, we consider the effect of different back- 
ground shapes in the fit. In Sect. 13 we adopt the 




r,„, (MeV/c^) 



Figure 4. The variation of the upper limit of 
Fee at 90% C.L. with respect to Ff. The slope of 
the solid line denotes the averaged ratio Fee/Ft = 
4.6 X 10"^. 



first order polynomial to depict the background, 
nevertheless, higher order polynomial can be used 
to delineate the background as well. 

In addition, we notice one background shape 
once adopted by DASP group, who tried to take 
into account the threshold effect of the charmed 
mesons |20| . 

6 

(Jdasv = CTe+e-^M+M- ' (^0 + ^ AjP^F^), (18) 

with 



l-s/(3.1 GeV )2 ' ^ ^ 

where describes the contribution from contin- 
uum; j ranging from 1 to 6 indicates the DD, 
DD*, D*D*, D,D,, D,Dl, D*D; thresholds, re- 
spectively. Aj are free parameters^; l3j are the 
velocities of the relevant particles; and F is a over- 
simplified form factor defined above. 

The fit results show that the effect of poly- 
nomial background is at the same level with or 

^It should be noticed that since the D*D* threshold has 
been taken into account in description of the Y (4260), the 
parameter Ag is set to be zero in the corresponding fit. 
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smaller than that of the DASP background, so as 
an estimation, we adopt the linear and the DASP 
backgrounds as two typical cases for background 
description. Comparing with those of linear back- 
ground fit, the results of DASP background fit 
for three summation schemes of -0 resonances are 
Tee < 50, 30, 250 eV/c2 at 90% C.L. respectively, 

Last, we also consider the possible effect of the 
fraction / on the measured Fee of the y(4260). 
Our fits indicate that with the increasing of factor 
/ the upper limit of Fee decreases and vice versa. 
The upper limit for Fee varies from 580 to 280 
eV/c^ when / changes from 0.3 to 0.7. 

Taking all the aforementioned possibilities into 
account, we adopt the most conservative result as 
our final estimation, that is 

pF(4260) ^ ^ (20) 

at 90% C.L. 

5. Discussion 

According to our study, the conservative esti- 
mation for the upper limit of Fee of the y(4260) 
is about 580 eV/c^ at 90% C.L., which is almost 
two times larger than the estimation presented in 
Eq. ©. 

Utilizing our upper limit Fee < 580 eV/c^ for 
the y(4260), together with the relation of Eq. (jTJ, 
we obtain the lower limit of the branching fraction 
at 90% C.L. to be 

S(y(4260) ^ TT+TT-J/i;) > 0.58% . 

This means that the partial width r(y(4260) 
TT+TT-J/ip) > 508 keV/c^ at 90% C.L., which 
is much larger than the corresponding par- 
tial widths of V' (89.1 keV/c^) [HI and Tp" 
(44.6 keV/c2) [aEa. 

CLEOc measured the cross section for 
TT+TT^ J/-0 channel to be 58 pb~^ jJS], which is 
consistent with the BaBar result, 50 pb""'^ Q]. As 
an estimation, we regard the central value cal- 
culated in Eq. to be the same for CLEOc 
and BaBar, but adopt the improved accuracy pro- 
vided by CLEOc, then we can obtain the follow- 
ing lower limits at 90% C.L. 

B{Y{A260) ^ TT+n-J/iP) > 0.6%, 



S(r(4260) ^ ttW/i/-) > 0.2% , 
S(r(4260) ^ K+R-J/iP) > 0.1% , 



and 

6(r(4260) ^ XJ/ii^) > 1.3% , 

where X denotes 7r+7r-(37), 7r07rO(8), K+K-{3), 
r?(5), 7r"(l), 77'(0), 7r+7r-7r°(0), and rir]{l). Here 
the numbers of the observed events by CLEOc 
are presented in parentheses. 

As we notice up to now only results on 
hidden-charm final states were reported about the 
y(4260); measurements involving open-charm are 
anxiously awaited to confirm existing specula- 
tions or provide clues for further theoretical in- 
quiry. In addition, more accurate Fee is still 
needed for a better understanding of the nature 
of the r(4260). 
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